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Abstract 



Recent results on jet production in neutral current deep inelastic scattering and high 
energy photoproduction at the HERA electron-proton-collider are briefly reviewed. The 
results are compared to QCD expectations in NLO and a s determinations using these 
data are summarized. 



5^ ; 1 Introduction 

o . 

The hard partonic interactions at the center of photon-induced hadronic processes 
in deep inelastic electron/positron-proton-scattering are charaterised by jets. The 
virtuality Q 2 of the exchanged photon can be varied from high values in deep inelastic 
scattering (DIS) to zero in photoproduction. According to the factorization theorem 
V ' the hadronic process can be factorized into a pertubatively treated hard part, which 

includes the evolution of the parton entering the hard process, and a soft part, 
which describes the parton distribution of the target nucleon; the characteristic 
energy for this separation is given by the factorization scale. In DIS processes 
involving light quarks only, an energy scale is provided by Q 2 or the transverse 



energy of the emerging hadronic jet (or jet system), whereas in photoproduction the 
jet energy is the only available energy scale; in heavy quark production the quark 
mass may provide an additional energy scale. In jet studies details of the hard 
QCD calculations and their matching to the evolution of the target partons are 
probed; the investigations may also provide information on the parton distribution 
functions (pdfs) of the proton and photon and can serve to measure a s and its scale 
dependence. 

In this report, recent results from the HI- and ZEUS-collaboration on the pro- 
duction of standard jets i.e. jets without heavy quark tagging in neutral current DIS 
processes and photoproduction are discussed; jets from diffractive processes and jet 
shapes are not included. 

*Talk given at the 9th Adriatic Meeting "Particle Physics and the Universe" in 
Dubrovnik/Croatia, 4.-24. September 2003 
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Figure 1: Leading order diagrams for jet production in photon-proton-interaction: (a) 
photon-gluon-fusion; (b) QCD-Compton process; (c) parton evolution 

For jet studies in DIS the Breit frame is preferred, which in the quark parton 
model (QPM) is defined as the reference frame in which the struck parton emerges 
collinear with the incoming purely time-like virtual photon q = (0, 0, —Q) such 
that 2xp + q = 0. Working in the Breit frame suppresses the QPM background 
and provides maximum separation between fragments from hard scattering QCD 
processes and beam remnants. Since the Breit frame is related to the photon-hadron- 
cms by a boost in proton beam direction (defined as +z-direction) , jet variables are 
convenient, which are boost-invariant - like the transverse energy Ex, the azimuthal 
angle <fi and the pseudorapidity rj = — hitan(9/2). 

The jet finding algorithm has to be stable against infrared radiation and collinear 
splitting and has to take into account the fact that at HERA not all final state 
particles are detected. In recent years, the algorithm preferred in jet studies by 
the HERA-collaborations has become the inclusive fcx-cluster algorithm pQ. The 
algorithms proceeds by calculating for each track or energy cluster i the quantity di = 
{Eti) 2 and for each pair i, j the quantity = mm(E Ti , E T j)[(r]i — f]j) 2 + (fa — 4>j) 2 ]- 
If of all resulting {0^,0^} the smallest one belongs to {di} it is kept as a jet and 
not treated further, if it belongs to {c%} the corresponding tracks or clusters are 
combined to a single object; the procedure is repeated until all tracks or clusters are 
assigned to jets. The algorithm is suited to handle overlapping jets, facilitates the 
separation between beam remnants and hard scattering fragments and yields smaller 
higher order QCD and fragmentation corrections than e.g. the cone algorithm. 

Jet cross sections are calculated by convoluting the matrix element of the specific 
hard scattering process with the pdfs of the target proton. Examples of LO hard 
scattering processes are boson-gluon-fusion and QCD-Compton scattering (fig. QJ. 
The parton evolution requires the summation of In Q 2 - and In 1/ x-terms, which is 
equivalent to a summing of ladder graphs pictured in fig.^ in the DGLAP-approach 
j2] only leading powers of lnQ 2 -terms are summed implying a strongly ordered 
increase of the virtuality of the exchanged partons towards the hard scattering. 
At small x where at HERA Q 2 is constrained to not too large values, the BFKL 
treatment of the parton evolution [3] is more appropriate, which resums In 1 /ar- 
teritis to all orders but retains the full Q 2 dependence; the parton chain is strongly 
ordered in x i.e. x -C x g <C ... -C 1, implying fc T -unintegrated (gluon) pdfs and 
off-shell matrix elements for the hard scattering process. In the CCFM ansatz [3] 
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colour coherence in the evolution ladder is taken into account imposing an angular- 
ordering of the gluon emission with the maximum angle determined by the hard 
scattering subprocess; the cross section factorizes into pdfs unintegrated both in k T 
and the gluon angle and an off-shell matrix element. At asymptotic energies CCFM 
evolution approaches BFKL at small x and DGLAP at large x and Q 2 . 

While at high Q 2 the virtual photon behaves as a pointlike particle, at low Q 2 
and in particular in photoproduction with Q 2 = 0, its hadronic structure has to 
be taken into account. The momentum transfer from the scattered electron to the 
parton cascade then takes place via a parton of the "resolved photon", which itself 
may evolve in a parton cascade towards the hard scattering process. Such resolved 
processes contribute a background to hard scattering processes but also can serve 
as a tool to investigate the photon structure function. 

The response of the detectors to jet events is studied with Monte Carlo gener- 
ated events, which are based on the measured pdfs of the proton and photon and on 
simulations of the fragmentation and hadronization of the partons. In RAPGAP the 
QCD matrix element is calculated in LO with NLO corrections evaluated by pertur- 
bative generation of parton showers; the hadronization of the partons is performed 
via JETSET, which is based on the Lund-string model. In ARIADNE, the parton 
cascade is generated using the colour-dipole-model (CDM) while in CASCADE it is 
calculated via CCFM evolution. 

2 Jets in DIS 

The inclusive jet cross section aj results from the convolution of the measured pdfs of 
the proton f a (x a , Hf) with the pertubatively calculated hard partonic cross section 
da a (x a ,a s ,fiR,fi F ), 

&j = J2 dxf a (x,/i F )da a (x,a s ,/i R ,/i F )(l + 8 hadl ) 

a J 

where the sum is to be taken over all partons a of the proton; aj has to be folded 
with the bremsstrahlungsspectrum. llr is the renormalization scale and lip the 
factorization scale. (1 + 5hadr) represents the nonperturbative correction accounting 
for the hadronization of the partons; it is estimated from MC generators. 

(i) Single jet cross sections have been measured by HI and ZEUS in the Breit 
frame using the fcx-cluster algorithm. The kinematical parameters for the HI data 
are 150 < Q 2 < 5000 GeV 2 , 7 < E T < 50 GeV, -1 < r? lab < 2.5, while for the 
ZEUS data they are Q 2 > 125 GeV 2 , E T > 8 GeV, -2 < ^ ab < 1.8. The good 
agreement of both measurements with NLO QCD calculations (DISENT) at high Q 2 
values (fig. EJ has enabled precision determinations of a s and its scale dependence 
from these data (see Sect. EJ). 

The ZEUS data have been used to study the azimuthal distribution of jets in 
the Breit frame [7j. The measured (^-distribution agrees well with the NLO QCD 
predictions from DISENT with either lir = or Q, providing an independent test 
of QCD. 
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Figure 2: Inclusive jet cross section in DIS in the Breit system (fig. from 6 ) 
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Figure 3: Inclusive jet cross section in DIS in the Breit system integrated over 5 < Q 2 < 
100 GeV 2 and 0.2 < r?i ab < 0.6, NLO QCD: DISENT without hadronization correction; 
shaded band: change of renorm. scale by 1/2 and 2 (fig. from |18j ) 
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H1 Forward Jet Data 
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Figure 4: Inclusive jet cross section in DIS in the forward direction i.e. 7° < Ojiab < 20°; 
CDM: Color-Dipole-Model; RG: RAPGAP; CASCADE: CCFM; (prelim, data, from 0) 



HI has extended the inclusive jet measurements [Hj to low Q 2 values of 5 < 
Q 2 < 100 GeV 2 , -1 < r/ lab < 2.8, for E T > 5 GeV and 0.2 < y < 0.7. The 
good agreement of the differential cross section in Et with NLO QCD (setting 

= Et) in the backward and central region worsens towards the forward direction 
i-e. f/iab > 1-5 and for Et < 20 GeV (fig. El)- in this region, NLO corrections are large 
and theoretical uncertainties due to the uncertainty in the renormalization scale are 
larger than the experimental errors. Further studies have shown the discrepancies 
to NLO QCD to develop at Q 2 < 20 GeV 2 . 

(ii) Data for inclusive jet production near the forward direction i.e. 7° < 8j\ a b < 
20° have been taken by HI [H] for 0.5 < Q 2 /E T < 2. The jet analysis using the 
&/r-cluster algorithm has been performed in the lab frame with cuts corresponding 
to 5 < Q 2 < 75 GeV 2 and 0.1 < y < 0.7. In fig. 0] the data are compared with 
different QCD models. Reasonable agreement with DGLAP is achieved only if 
resolved photon processes are included, but also the colour dipole model describes 
the data reasonably well. 

A process related to forward jet production is forward 7r°-production, which has 
been studied by HI [TU]; here the forward parton is tagged by a single energetic 
fragmentation product i.e. the 7r°. While in principle smaller angles than in jet 
production can be reached, the cross sections are lower and the hadronization un- 
certainties are higher. The data can be described by DGLAP models only after 
including resolved photon processes, but also BFKL models provide a reasonable 
description of the data. 

(iii) LO contributions to dijet production are the QCD-Compton effect (QCDC) 
and boson gluon fusion (BGF) (see fig. Q). In the high Q 2 region where QCDC is 
dominating and the pdfs are well constrained by fits to F 2 data, dijet measurements 
using ot s as input can serve to test pQCD; alternatively, by comparing the data with 
pQCD, a s (and its running) can be determined. At low Q 2 , where the cross section 
is dominated by BGF, the comparison of data with assumptions on the pdfs can 
constrain the gluon density distribution. 
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Figure 5: Inclusive jet {da tot) and dijet (afo"2+i) cross sections with dijet fraction R2+1 in 
DIS; shaded band: uncertainty in absolute energy scale of jets (fig. from [TT] ') 



Dijet data have been taken by HI pQ and ZEUS ^JJ in similar kinematical regions 
and analyzed by similar methods. Asymmetric cuts have been applied in Et in order 
to avoid infrared sensitive regions of the phase space. At not too low Q 2 and Et 
i.e. Q 2 > 10 GeV 2 and E T > 5 GeV, the data are reproduced (fig. |3J) within 10% 
by NLO QCD (DISENT). For the a s analysis of these data see Sect. El 

Recent dijet-data from HI in the low- a; region i.e. 10~ 4 < x < 10 -2 for 
10 < Q 2 < 100 GeV 2 agree within errors with DGLAP based NLO QCD predictions 
for inclusive cross sections. Serious disagreement is observed, however, in the ratio 
S of the number of dijet events with azimuthal separation <fi < 2n/3 to the total 
number of dijet events. While the LO (DISENT) prediction is much to low, the 
agreement is improved except in the low-x and low-Q 2 region if a third hard jet is 
included in the NLO (NLO JET) calculation (fig.®. A comparison with RAPGAP 
shows that the inclusion of resolved processes improves the agreement except in the 
very low-x and low-Q 2 region. From a comparison of the data with the CCFM 
based CASCADE MC a strong sensitivity to the assumptions on the unintegrated 
gluon distribution can be inferred, which may be used to constrain the gluon pdf. 
ARIADNE, which is based on CDM, describes the low- a; and -Q 2 region reasonably 
well, fails, however, at larger Q 2 values. 

The jet studies in deep inelastic scattering have been extended to trijets [T3j . 
[Tlj as well as to subjets [To] . 

3 Jets in photoproduction 

The description of photoproduction involves the convolution of the hard partonic 
scattering cross section with the pdfs of both the proton and the photon. In standard 
jet production i.e. jets from light quarks only, for the energy scale hr = hf = H = 
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Figure 6: Ratio S of number of dijet events in DIS with azimuthal sep. 4> < 2n/3 between 
the two highest transverse momentum jets to the total number of dijet events compared 
with LO and NLO QCD predictions for dijet and trijet production; light shaded band: 
quadr. sum of hadroniz. error (dark shaded band) and renorm. error; (fig. from |12j ^ 



Et is taken. Apart from the folding with the electron bremsstrahlung spectrum, the 
jet cross section is written as 

Uj ^ ' J J dx~fdXpfp a (Xp, /i ^fjbfaji A* ^d&ab^Xp, X*y, Oi s , ) ( 1 -|- 5hadr ) > 
a,b 

where the sum is to be taken over all partons a, b from the proton and photon. 

(i) In single jet production a wider kinematic range is accessible than in dijet 
production, cross sections are higher and systematic errors due to the treatment of 
infrared unsafe regions are avoided; on the other hand, the reaction is less sensitive 
to details of the hard scattering process. 

In recent measurements at Q 2 < 1 GeV 2 of differential cross sections by ZEUS 
[TH] and HI using the /cx-cluster algorithm in the lab frame, the kinematical 
range of E T > 5 GeV, — 1 < r]i ah < 2.5 and the photon-proton-cms energy W of 
95 < W < 293 GeV has been covered. While LO QCD calculations fail to reproduce 
the data in their normalization, good agreement is found with the corresponding 
NLO predictions (fig.Q left). A subset of data taken by HI at Q 2 < 0.01 GeV 2 for 
5 < Et < 12 GeV indicate an excess over predictions towards increasing 7] values. 
The precision of the experimental data and of the predictions is not yet sufficient 
to discriminate between different proton and photon pdfs. For the determination of 
a s (M z ) and its scale dependence from the data of [T2] see Sect. El 

In the ZEUS experiment ^Hj also the scaled invariant jet cross section 
E T Ed 3 a /dp x dpydp z has been measured as a function of the dimensionless variable 
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Figure 7: Inclusive jet production cross section in photoproduction. Ratio of scaled 
invariant cross sections (av. over —2 < < 0) for the two shown (W 7P ) compared with 
NLO QCD; (fig. from [E]) 

x T = 2Et/W 1p at (W 1P ) = 180 and 255 GeV. The cross section ratio as function of 
xt violates scaling as predicted from QCD (fig. right). 

(ii) The dijet photoproduction cross section is dependent on the dijet an- 
gle 9* in the parton-parton-cms, which is sensitive to the matrix element of the 
hard scattering subprocess. Direct photon processes are dominated by BGF with 
the angular dependence of the spin 1/2-progator ~ (1 — cos#*) -1 whereas in re- 
solved processes gluon-exchange with the steeper angular dependence of the spin 
1-progator ~ (1 — cos#*)~ 2 is dominating. The dijet events can be enriched with 
direct or resolved processes by cutting on the fraction x° bs of the photon-momentum 
transferred to the two jets of highest transverse energy, which is given by x° bs = 
{Erpie~' ni + Et2G~ V2 ) /2yE e with y the fractional electron energy carried by the pho- 
ton in the p-rest frame. 

Dijet photoproduction has been studied by ZEUS [IB] and HI [THj under similar 
kinematical conditions using the A/p-cluster algorithm. The differential cross sections 
in the dijet invariant mass, and 77 agree well with NLO QCD expectations and 
in particular the measured dijet angular distributions (fig. |BJ) shows the expected 
difference for direct (x° bs > .75) and resolved processes (x° bs < .75). The sensitivity 
of the data to the parametrization of the photon-pdfs appears to be affected by the 
choice for the £ T -cuts (ZEUS: E T1 > 14 GeV, E T2 > 11 GeV; HI: E T1 > 25 GeV, 
E T2 > 15 GeV). 

4 ^-determinations 

The good agreement of jet cross sections with QCD predictions has enabled deter- 
minations of a s and its scale dependence. 

The HI inclusive jet cross sections jS] as a function of have been fitted to 
the QCD prediction for four regions of Q 2 using the proton-pdfs and /jlr and /jf as 
input. The resulting a s has been shown to be stable against variations of the jet 
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Figure 8: Inclusive dijet photoproduction as a function of the dijet angle 0* in the parton- 
parton-cms for different cuts on the fractional photon momentum x°^ s transferred to the 
dijet system. Comparison with NLO QCD for different photon-pdfs; (fig. from 18 ) 



algorithm. The combined fit to all Q 2 -data evolved to mz is shown in the summary 
fig. El). 

For the analysis of the ZEUS data [B] the following procedure has been applied: 
The cross sections da/dA with A = Q 2 , Et have been calculated in NLO QCD for 
the same pdf-set with three a s -values. These calculations were used to parametrize 
the a s -dependence of the binned cross sections (da/dA)i for each bin i according to 
(da(a s )/dA)i = Cu. a\ + C 2 i- a 2 with a s = a s (M z ). From a x 2-n t of this ansatz 
to the measured da/dA, a s was obtained for the chosen regions of A. The best fit 
was obtained for Q 2 > 500 GeV 2 (fig. EJ). The same procedure has been applied to 
determine the scale dependence of a s . With Et as energy scale, the a s dependence 
of da/dEx was parametrized in terms of ol s {(Et)) where (Et) is the mean value of 
Et in bin %. The result is shown in fig. 

The same method has been applied to other jet results from ZEUS i.e. the 
dijet fraction jTJ], subjet multiplicities (jTHj) and jets from photoproduction [TH]. A 
summary of the results is shown in figs. IHlandlTUl 



5 Summary 

In the analysis of jet production in DIS and photoproduction at HERA the lon- 
gitudinally invariant A>r-cluster algorithm in its inclusive version has become the 
standard jet finding algorithm. 

In DIS, the inclusive jet cross sections at higher values of Q 2 and Et are well 
described by NLO QCD; in this region higher order and hadronization corrections 
are small; this holds as well for photoproduction. 
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Figure 9: a s (Mz) from recent HI and ZEUS measurements and the world average 
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Figure 10: Scale dependence of a s from recent HI and ZEUS measurements and QCD 
predictions for a s (Mz) = 0.1183 (world average) 
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In the forward region and/or at smaller Q 2 resp. x values, the situation is less 
satisfactory. DGLAP based calculations are expected to become less reliable in 
this region, the mentioned corrections are sizable and the hadronic structure of the 
photon i.e. resolved processes have to be taken into account. Calculations based 
on BFKL or CCFM evolution only partly show better agreement. An increase in 
experimental statistics and efforts to reduce the theoretical uncertainties are highly 
desirable for a better understanding of this challenging region. 

The measurements of a s and its running from jet data have yielded results which 
are of perhaps unexpected precision; they can well compete with other precision 
measurements e.g. from e + e~-annihilation and are in good agreement with the 
world average. 

Acknowledgements I thank G. Grindhammer for valuable discussions and com- 
ments on the paper. 
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